Universal "rephasing" dynamics 
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We consider a quantum quench in which two initially independent condensates are suddenly coupled, and 
study the subsequent "rephasing" dynamics. For weak couplings the time-evolution of physical observables is 
predicted to follow universal scaling laws, connecting the short-time dynamics to the long-time non-perturbative 
regime. We first describe a two-modes model valid in two and three dimensions and then move to one dimension, 
where the problem is described by a massive Sine-Gordon theory. Combining analytical and numerical meth- 
ods we compute universal time-dependent expectation values and distribution functions, allowing a qualitative 
comparison with future experiments. 



Equilibrium systems near second-order phase transitions 
are the best known examples of universal behavior in many- 
body systems. For these systems, universality implies that all 
macroscopic quantities are related to the distance from the 
transition point through universal critical exponents. In re- 
cent years, many efforts have been dedicated to the extension 
of these concepts to quantum systems out of thermal equilib- 
rium both in openflJQ and closed systems. In the latter (see 
Ref. [4] for a recent review) some authors considered systems 
that are slowly driven through a second-order transition and, 
extending the Kibble-Zureck argument, found universal scal- 
ing laws of static and dynamic quantities||5}{8). Others con- 
sidered sudden quenches and predicted a universal scaling of 
physical observables in the steady state regimell9 UTD . Here 
we extend these arguments by demonstrating examples of uni- 
versal behavior in the transient dynamics following a sudden 
quench. 

Here we consider systems that are initially prepared in the 
ground state of a gapless Hamiltonian. The quench corre- 
sponds to the sudden opening of an excitation gap A (see 
Fig[TJa)). In the Renormalization Group language, this cor- 
responds to the sudden switching-on of a "relevant" pertur- 
bation. The resulting dynamics is described by a strongly 
interacting theory and is intrinsically non-perturbative. Ex- 
perimentally, this can be realized for example by preparing 
two independent Bose-Einstein condensates in their respec- 
tive ground state and then suddenly turning on a finite tunnel- 
ing coupling j x (cf. Fig[TJb)). This opens an excitation gap 
for the relative-phase modes, given, at a mean field level, by 
the plasma frequency A = 2 ^Jfij x , where fi is the chemical 
potential. Beyond-mean-field effects will be extensively dis- 
cussed below. 

In the limit of j x <§c ft, the dynamics following the quench 
is dominated by the low-frequency modes and is thus expected 
to be universal. For a generic physical observable C(f) we 
predict the universal contributions to have the form 

C(t) = ij) R(At). (1) 

Here rj is the scaling dimension corresponding to the spe- 
cific observable and R is a scaling function, fixing the uni- 
versal behavior of the transient many-body dynamics. As we 
will demonstrate, the scaling ansatz ([T]i is valid at times both 
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FIG. 1: [Color Online] (a) Dispersion relation before (dashed) and af- 
ter (solid) a quantum quench opening a gap A and involving only the 
low-frequency modes, (b) Physical realization: sudden coupling of 
two independent condensates via a uniform tunneling j ± . (c) Lattice 
model used for the numerical calculations. 

shorter and longer then the inverse gap 1/A. This allows to 
extract information about non perturbative steady states from 
the short-time dynamics, where time-dependent perturbation 
theory applies. A similar relation is known to hold for clas- 
sical systems, where it has found important theoretical and 
computational consequences IT2"lfl5l . 

Two modes model - To substantiate our predictions we 
first apply a two-modes theory, valid in two and three 
dimensions [49 1 . By considering only the macroscopically 
occupied state of each condensate, we describe the sudden 
quench byEj) 

H(t) = £ J] 6n « + j Jt) W^ 2 + Hc ) ■ (2) 

Here \p\ and fa are canonical bosonic operators, 5n a = ij/ a ij/ a — 
N a , and N a — Np — N/2 is the average number of particles 
in each condensate. The two terms in Eq.Q describes re- 
spectively the interactions among the atoms in the same con- 
densate and the sudden switching-on of a tunneling coupling: 
j±(t) = j±®(t), where 0(f) is the Heaviside step function. We 
postpone the discussion on the conditions of validity of this 
model to the end of the paper. Exact diagonalization of Eq.Q 
is shown in Fig. |2ja), for different values of j x . Subplot (b) 
shows the excellent data collapse obtained by rescaling the 
time axis, and demonstrating the validity of our scaling ansatz. 
To describe the observed scaling behavior, it is useful to 
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FIG. 2: [Color Online]Time evolution of the interference contrast 
C 12 (0 = <^JiA2 + H.c>/iV = 2cos(V20), according to the two- 
modes quantum Hamiltonian |2]( with N = 1000. (a) Exact diago- 
nalization for different values of the ratio j±/fi = 0.10 (red), 0.060 
(green),0.037 (blue), 0.022 (cyan), 0.0135 (magenta), (b) Same plots 
on rescaled axes, showing the universal behavior ^ with rj = 0, su- 
perimposed to the analytic solution of the simple pendulum (thick 
solid line). 



introduce the phase fields <p a , canonically conjugate to 5n a , 
through the approximate relation iff a * VA^e* . Moving to 
the relative coordinates <p - {<p\ — 4>2)l V2 and n = {5n\ - 
5ni)j y2, we than obtain the Josephson junction Hamiltonian 



H(t) 



2j" 2 
— n 
N 



20(f) j ± N cos( V20). 



(3) 



The low-energy states of this Hamiltonian are well described 
by an harmonic approximation, in which the cosine is ex- 
panded to quadratic order in </>. For the sudden quench, this ap- 
proach would predict endless oscillations. Instead, as shown 
in Fig(2] the universal dynamics displays strongly damped os- 
cillations, with a damping coefficient of order one. The origin 
of this large damping is due to the higher-order Taylor compo- 
nents of the cosine. In the initial state, | if/(t — 0) ) = | n = ) , 
the inital phase is completely undefined and all terms con- 
tribute to the dynamics by the same amount, highlighting the 
strongly interacting nature of the predicted dynamics. 

In the limit of large N, an alternative, non-perturbative 
description is available: the semiclassical approach, or 
Truncated Wigner approximation I! 16ifT8ll . Following this 
method [51], we map the quantum Hamiltonian pi to the 
equations of motion of a classical simple pendulum cr<f>/dt 2 = 
-A 2 sin( y2(f>)/ V2, whose analytical solution is known. The 
quantum nature of the problem enters through the initial con- 
ditions, for which </> is uniformly distributed between zero and 
V2~7r. As shown in Figjijthis approach describes damped os- 
cillations with a single frequency scale A and correctly repro- 
duces the numerical solution of the full quantum model. 

One dimension - In one dimension the initial state is not 
described by two macroscopically occupied modes, but rather 
by two "quasi-condensates" with power-law correlations. Us- 
ing the standard bosonization technique] 19 1, we describe the 
system in terms of two phase fields, (f>i(x, t) and <pi(x, t), and 
the time-dependent Hamiltonian 



H - ^ H LL [<p a ] - 2jjt)p | dx cos(^ I - 



<h)- (4) 



two independent condensates as "Luttinger liquids" ll20ll with 
sound velocity v and Luttinger parameter K. The latter mea- 
sures the ratio between the longitudinal kinetic energy and 
the interaction energy: K — > oo corresponds to free bosons 
and K — 1 to the Tonks-Girardeau limit. The second term of 

Eq. (HI describes the tunneling between the two condensates 
— t t 

JjAOwi^+i^W 1 )' as obtained through the "bosonization dic- 
tionary" 01,2 — > y/poe"^ 1 - 2 . 

Introducing the relative coordinate <p = (4>\ -02)/ V2 we 
than obtain the well-known "Sine-Gordon" model with a non- 
linear term proportional to cos(0! - fo) = cos( V20). This 
model is characterized by a quantum phase transition at the 
universal value of K = 1 /4, separating a gapless phase (K < 
1 /4) from a gapped phase (K > 1/4). 

The effects of the quench in the two phases are completely 
different. For K < 1/4, the cosine is irrelevant and excites 
primarily the high-frequency modes ||2T1|221 . Accordingly, it 
is natural to expect the resulting dynamics to be non universal. 
This situation is reversed for K > 1/4, where the cosine is rel- 
evant and opens a gap[52|. In this case the average amount of 
energy per mode scales to zero much faster than the gap ||23l 
and, as we will now show, the dynamics of the system is dom- 
inated by quantum effects. 

To study the dynamics of the problem we perform a non- 
perturbative expansion in the coupling strength j x . We first 
move to a path-integral formulation of the problem (For an 
introduction see for example Ref. EH [25]) and represent the 
Hamiltonian Q is terms of the real-time action: 

S = f dt fdx- [(d x (f> 2 ) - (d t <f> 2 )] + 20(f) j ±Po cos( V2>) 

(5) 

Here jk is the Keldysh contour and we switched to units 
where v = 1, Splitting the Keldysh contour in its forward and 
backward branches (<f>+ and 0_ respectively) and then moving 
to their symmetric and anti-symmetric combinations (defined 
by <p ± = (</>+0)/2 and often termed "classical" and "quantum" 
components), we have: 

2 I dt cos( V20) = 2 I dt cos( V2</> + ) - cos( V20_) 

Jy K J 

= 4 J dt sin( V2</>)sin( V20) = ^ J dt 



e=±, l=± 



(6) 



Expanding e' s in a Taylor series around j ± = we obtain 
Cn(t) = <cos( V2>(f) )> 



(ijxPof 



dx e 



V2;'(e0+e 



> 

(7) 



Here the first term describes the low-frequency modes of the 



Note that the A^-th term of |7| contains 2N continuous inte- 
grals and 2N discrete sums, whose dummy indices we will 
denote as (x,, f,, e,, e,) with i = l,...,N. In analogy to the 
equilibrium case l26ll27ll . we can interpret the first three ele- 
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FIG. 3: [Color Online]Interference contrast obtained from the series 
expansion |9](. Solid curve and error-bars: numerical estimation of 
the universal curve, obtained by computing N m!a terms of the ex- 
pansion using 10 6 random evaluation points. Dashed curve: analytic 
result for the first order contribution N = 1, obtained in the S. M.. 
Numerical values: (a) K = 1.2 and iV max = 11; (b) K = 5 and 



(b) 




FIG. 4: [Color Online] Interference contrast obtained from the Trun- 
cated Wigner approach (TWA) for K = 5 (System size: L = 400, and 
number of random paths:,/V = 400). (a) Time evolution of (e 2 *") for 
different values of jjyc= 0.27 (red), 0.16 (green),0.10 (blue), 0.060 
(cyan), 0.036 (magenta), (b) Data-collapse according to (10) , super- 
imposed to the results of the series expansion (black curve, identical 
to Fig0») 



(a)K= 1.2 



(b) K = 5 



(a) 



merits as the space-time coordinates and charges of a Coulomb 
gas with N particles (dual to the atoms of the original prob- 
lem). In contrast, the e, indices do not have an equilibrium 
analogue and rather impose a light-cone[53|, constraint limit- 
ing all contributing "charges" to the past-light-cone of (x, t), 
or |x, - x\ <t - tj. The relevance of the "light-cone physics" is 
generic to global quantum quenches (see for example Refl28l 
and prevents an exact mapping to an equilibrium effective 
model. 

The expectation value appearing in Eq. (|7]l refers to the 
initial state, corresponding to the ground state of two inde- 
pendent Luttinger liquids, and can be computed using the 
quadratic part of Eq. Q. After rescaling the time and space 
variables of integration, we arrive to 



Cl 2 (f) = 2 CW ( ^ f) 



\2N-(N+V)/(2K) 



U(AK-l) 



2K/(4K-l) f \2N-(N + l)/{2K) 



(8) 



(9) 



The coefficients are universal! 54 1: they do not depend nei- 
ther on the cutoff nor on j ± . 

To achieve a universal scaling of the form of Eq. ([TJ by 
rescaling the time in units of the inverse gap, A -1 . For the 
Sine-Gordon model, the excitation gap was obtained using the 
form-factor approach in Ref 29jand, in our notations, it is pro- 
portional A = fi(j_i_/p) 2K/{4K ~ l \ After this rescaling, Eq. ^ 
becomes 



R(At) (10) 



with rj = l/(2K). As expected, in the limit of non-interacting 
bosons K — > oo we restore the mean field result presented 
above. Remarkably, the same anomalous scaling dimension rj 
obtains in the ground state of the final Hamiltonian[30], high- 
lighting the quantum nature of the present problem. 

The above calculations were obtained by using the asymp- 
totic values of the correlation functions and thus neglect the 



non-universal contributions due to the high-frequency modes 
of the system. At times larger than the inverse cutoff t » l/p, 
these modes can contribute a non-universal constant factor, 
which may hide our universal predictions. Being related to 
the high-energy physics, the non-universal contribution must 
be an analytic function of j ± and can, at most, be a linear 
function of j±/p. As we saw, the universal contribution scales 
as f ± it K ~V anc [ j s therefore dominant for any K > 1/2. Re- 
markably, the marginal case of K — 1/2 can be exactly solved 
by mapping the problem to non-interacting fermions PD . At 
this point, the universal and non-universal contributions scale 
in the same way and the latter are no more negligible. 

At present, we were not able to compute the universal func- 
tion R by the analytical resummation of Eq. i9). Instead, we 
can numerically estimate it by averaging over random sam- 
plings of (e,, Xu f,). The resulting universal curves are shown 
in Figj3] for two different values of K. As expected, this ap- 
proach is limited to time scales of order 1/A due to a "dy- 
namical sign problem", here related to the alternating signs of 
the coefficients c^. To predict the universal behavior at longer 
times, we will now present two alternative numerical methods. 

The first method extends the semiclassic approach pre- 
sented above for the two-modes model. The corresponding 
classical equations of motion are: 



3 2 = 3 2 0-^;> o sin(V2>). 



(11) 



Fig ^demonstrates a good agreement between the numerical 
solution of Eq.( 1 1 1 and the above non-perturbative expansion. 
Note however that Eq.([TT]i is characterized by a single energy 
scale y/2nj ± po/K ~ jY\ matching the excitation gap of the 
original quantum model only for K » 1. As expected, the 
semiclassical approach breaks down for K ~ 1, when the sys- 
tem approaches the quantum phase transition. 

To describe the strongly interacting regime K ~ 1, we map 
the problem into a lattice model and apply one of the deriva- 
tives of the Density Matrix Renormalization Group (DMRG) 
algorithm ll32l [331 . This approach is well suited to situations, 
as the one under present consideration, in which the time evo- 
lution is performed with a gapped Hamiltonian. To minimize 
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FIG. 5: [Color OnlineJInterference contrast obtained from the nu- 
merical solution of the hard-core bosons model \12\ with V/J = 
-1/2, corresponding to K = n/(2n - 2 arccos(V/2/)) « 1.2. (a) 
Time evolution of Cnif) for different values of j±/J= 0.27 (red), 0.16 
(green),0.10 (blue), 0.060 (red), 0.036 (cyan), 0.022 (magenta), (b) 
Data collapse according to ^10\ , superimposed to the results from the 
series expansion (black curve). 

Methods: The calculations where performed using the Time Evolv- 
ing Block Decimation (TEBD) algorithm of Ref.[36 37 1 for system 
size L = 100, number of states x = 100 and took around three 
days/curve on a single core. With respect to Fig|3[a), the black curve 
has been rescaled to account for the mapping from the microscopic 
model \\2\ to the effective field theory l|4j. 



the memory requirements we consider hard-core bosons on a 
ladder (see Fig{TJc)) represented by the Hamiltonian 

n r 



+ j±(t)(b{ i b 2 j+n.c.) 



(12) 



Here the energy scale j\\ sets the chemical potential and V < 
is a nearest neighbor attractive interaction, meant to partially 
counter-balance the hard-core constraint (55]. For j ± = the 
two independent chains are Bethe-ansatz[34 | solvable, allow- 
ing to exactly determine the Luttinger parameter! 20, 33. The 
results of this calculation for the expectation value Cnit) are 
presented in Fig(5] and confirm our scaling ansatz. (See the 
caption for technical details.) 

Experimental realization - In the above derivation, we con- 
sidered the ideal limit of a system initially prepared at zero 
temperature and suddenly quenched. In actual experiments, 
one necessarily deals with finite initial temperature 7q and a 
finite ramping time r. Intuitively one would expect the ideal 
limit to be valid as long as the temperature and the ramping 
time are smaller than the (inverse) frequency of the highest 
excited state reached by the quench. For the Hamiltonian 
Q, this energy scale corresponds to Nj ± = N exc i te dA, with 
Abated = N yj v'j_/yU » 1. In contrast, the numerical ll56l and 
the semi-classical analysis[57 1 show that the only relevant en- 
ergy scale is A, in agreement with our claims of universality. 
The solution of this apparent contraddition is related to the 
many-body nature of the problem: As discussed in Ref.|38|, 
the most probable excitation path does not directly connect 
the ground state to the highest excited state, but is rather given 
by subsequent steps connecting neighboring levels, whose en- 
ergy splitting is of order A. 

To determine the feasibility of the proposed experiment, 



one should combine the two above constraints, Tq <k A and 
t <k 1/A, with the weak coupling condition \!j±/fi — A In «: 
1. The data collapse observed in the right panels of Fig.s|2]|4] 
and[5]demonstrate that the universal scaling can be already ob- 
served for j ± < O.ljU, or A < 0.5/i. Temperatures (ramp times) 
smaller than the (inverse) chemical potential by one order of 
magnitude seem to be reasonable assumptions for present ex- 
periments with neutral atoms and would allow the abservation 
of the predicted universal dynamics. 

The most natural realization of the proposed experiment in- 
volves wave interferometers of atoms on a chip |t39l . This re- 
alization poses a serious challenge for the observation of one 
dimensional physics, due to the extremely high values of the 
Luttinger parameter K * 50. In this limit the average co- 
herence Ci2(f) is mostly dominated by mean field dynamics 
(compare for example Fig|2]with Fig|4|. One possible solu- 
tion to overcome this difficulty is to study the spatially inte- 
grated coherence|27j. In the mean field model Q this ob- 
ject is a pure phase with unit modulus. Quantum fluctuations 
along the longitudinal direction determine a broader distri- 
bution. At t = the two condensates are independent and 
the distribution function is known to correspond to the Gum- 
bel distribution l27l l40l ATI . During the time evolutionll58l. 
we observe that P(\Z\ 2 ) oscillates between the initial Gum- 
bel distribution and narrower peaks, indicating the forma- 
tion of a long-range order in the longitudinal direction. Fol- 
lowing the arguments given before and rescaling both time 
and space coordinates by 1 / A, it is possible to show that the 
full distribution function of Z follows the universal relation 
P(Z) = P(Z, At, AL), allowing to perform a combined finite- 
time and finite-size scaling analysis. 

It is instructive to compare our findings with the experimen- 
tally observed "dephasing" process in which one condensate 
is suddenly split in two parts Il27l [39l I42H441 . In this case, the 
time evolution is performed with a gapless Hamiltonian and 
the resulting dynamics can be studied by the exact solution 
of the quadratic Luttinger Liquid Hamiltonian (Eq. Q with 
j± = 0). Due to the absence of an excitation gap, at long-times 
the "dephasing" process is dominated by a thermal rather than 
quantum behavior] 44 45]. In contrast, in the present "rephas- 
ing" experiment we predict a universal dynamics dominated 
by quantum effects both at short and long times. 

In summary, we studied the many-body quantum dynam- 
ics of two initially independent condensates, suddenly cou- 
pled via uniform tunneling ("global quench"). In two and 
three dimensions we applied a "two-modes" model, while in 
one dimension we described the problem in terms of a Sine- 
Gordon theory. In both cases, the tunneling term is relevant 
and the observed dynamics is governed a non-perturbative 
strongly interacting theory. Neverheless, we find that the dy- 
namics at both short and long time scales (with respect to 
the inverse gap) is described by the same universal scaling 
laws. This may have important implications for both numer- 
ical calculations and experiments, in which long-time evolu- 
tion is hard to achieve. It may also allow to perform finite- 
size scaling and renormalization-group analysis in the time 
domain ll2Tl l46l l47l . in analogy to the established equilibrium 
momentum-shell techniques. 
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Appendix A: Semiclassical approach to the mean-field 
Hamiltonian (3) 



Consider the generic quantum time-dependent Hamiltonian 

(Al) 



H(f) = ^/+2JV; i /(0V(i), 



where x, p are unitless canonical variables, and V(x), /(f) 
are unitless functions. For the ramp of a Josephson junction , 
V(x) = cos(x) and 



fit) 



The correspondent action is: 



t/T t < T 
1 t > T 



(A2) 



-/ 

•JJk 



S [ .v. p 1 = S„ + i \ dt pd t x -^P 2 - 2Nj\fW(x), (A3) 



where jk runs over the contour (0 — > t — > 0) and e' s ° = 
( iffQ | pt,xt ) ( Pq , Xq \ ) sets the initial conditions. If the 
initial state is | p = ) , then we have 



e is ° = 6(p + )6( Po )6(x + - x-) 



(A4) 



After integrating out p you get: 

N 



S[x,p] =S + iN f dt^ 7 (d t x) 2 -2j ± J{i)V{.\) (A3) 

•JJK 



In the limit N — > oo, only the saddle point contributes to the 
path integral. As always, the saddle point, corresponds to the 
classical equation of motions 



d]x = 2co f(t)V'(x) 



(A6) 



with initial conditions substituted by the distribution function 
P(x ,f ) = e iSo[XoAxo] . Note that the energy scales J = ; X 7V 
and U — p/N have completely disappeared from the problem. 



Appendix B: Technical details of the nonperturbative expansion 

Eq. (8) 



The starting point of the present derivation is Eq.(7). If we 
define </>o = <p(x, t), eo = €q = 1 and use the compact notation 

= and iv = Uti IT dt > fZ, dx >> we can ex P ress 

the intereference contrast as 



Ci2(r)=£ 



2 N N\ 



£& <exp 



2i V 6i(pi 



\ i=0 



+ en 



> (Bl) 



The expectation value refers to the initial state, given by 
the ground state of two independent Luttinger liquids, and 
can be computed from the action (5) with j x = 0. Deal- 
ing with a quadratic action, we can apply Wick's theorem 
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(e iA ) = { e - (Al) l 2 ) and obtain 
n , A y (ihPof f , 



exp 



,2 \ 



) 

(B2) 



Yj lJ± N\ I £ ' 6XP ~ Tj 



V ij=0 



x e, exp 



(B3) 



where we used the identity ($> 0) = 0, related to the causality 
constraintE25ll. 

In one dimension, the correlation functions ((pi4>j) are infra- 
red divergent. As a consequence the above expression van- 
ishes, unless Yif=o € i - 0- ^ we interpret e, as the charges of 
a Coulomb gas expansion, this condition coincides with the 
requirement of a total charge neutrality. It constraints N + 1 
to be even and Yli e i - (-1) <A,+1 ^ 2 . This alternating sign will 
be the origin of a "sign problem", preventing the numerical 
resummation of the series. 

We now assume that the largest contribution to integrals 
over time and space comes from the universal asymptotic be- 
havior of the correlation functions [48 1 



(B4) 



(((Pi - (f>j) 2 ) = — log (yU - Xjf - (ti - tjf 

(<f>ih) = ^®( t i- t j-\xi- x j\)) 



We discussed the non-universal short-distance behavior at the 
end of the paper and show that they indeed are negligible for 
K> 1/2. 



Using ( |B4| >, we obtain that the second line of ( |B3| > becomes 

f n \ 



£j exp 



!=0 



= -2; sin 



(B5) 



Here we explicitly performed the sum over e, and we defined 
rij as the total charge of the points inside the future light cone 
of (xj, tj). Note that this expression vanishes unless nj + 
for all j = l,.,,N. This condition can be fulfilled only if all 
the "charges" are inside the past light-cone of (xq, to) = (0, t), 
satisfying | jc — jcy | < t—t;. This "light-cone physics" is relevant 
to any global quantum quench (see f.e. Refl28l and highlights 
the difference between similar equilibrium Coulomb expan- 
sions. 



The first line of (B3 i equals to 



exp 



2^^;*W = [] \(Xi-Xj) 2 -(ti-tj) 
K i,j=0 ' i,j=0Ji=j 



(B6) 



Note that there are (N + 1)N homo-signum combinations 
and (N + 1)(N +1) etero-signum combinations, thus at long 

. / 1/7Jr \(Ar+l)JV/2-(Ar+l)W+l)/2 

times this term is proportional to It ' 1 = 




T Q /u 



FIG. 6: Comparison between the idealized case considered in 
Fig.2and (a) finite initial temperature To, or (b) finite ramp time t. 
The solid lines represent 6C l2 (t) = Cff c (r) - Cj d 2 cal (f) at different 
times t = 10/yU (red), 20/fj (green), and 30/i (blue). The vertical 
dashed lines correspond to (a) Tq = A, (b) t = 1 /A. As predicted, 
the dynamics on the left side of these lines is independent of T and 
t. Numerical values: = 0.01, A = lyjJZJJi = 0.2/i, and 

N = 1000. 



t ( N + l )/( 2K \ To obtain Eq. (8) it is sufficient to rescale the time 
and space integration parameters dt — > t dt and dx — » t dx. 

The coefficients are universal (do not depend neither on 
the cutoff nor on the gap) and are given by 



s 



2 JV (-l) ! r 
AH 



N+l M 



dt; 



sin(-n,— ) Y\ ~ x jf ~ ft - tj) 



2|^,./(4X) 



(B7) 



with (xo,to) = (0,1) and eo = 1- The coefficient c\ corre- 
sponds to the first order time-dependent perturbation theory 
and can be computed analytically 



c\ — 2 sin 



dx 



1 



| x 2 _ t 2\l/(2K) 



2sin(^) V^r(l-^) 

( 2 -^) r (I-^) ' 



(B8) 



(B9) 



Here the second line has been computed with Wolfram's 
Mathematica (v8.0.1) and F is the Gamma function. 



Appendix C: Finite temperature and finite ramp time 

We now consider the effects of a finite initial temperature Tq 
and a finite ramping time r. Because the dynamics described 
in this paper is characterized by a single time scale A -1 it is 
natural to assume that the universal dynamics should not de- 
pend on these parameters as long as To <K A and r <K 1/A. 
For the mean-field model (3)this claim can be verified by ex- 
act diagonalization, as shown in Fig|6] or by the semi-classical 
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method introduced above. 

For the one dimensional case, one can move one step fur- 
ther and show that a finite Tq and r can be used as additional 
scaling parameters: 

C12W = <cos( V20(O)> = (xf*( A '' At ) (CI) 

In particular, this scaling relation confirms that, as long as 
Tq <k A and t «: 1/A, the resulting dynamics should be in- 
dependent on these parameters. To demonstrate the validity 



of Eq.lCl 1 and starting from the non-perturbative expansion 
presented above, one obtains: 

Ci2 ^ = Z ^SSt- f eiiiF (-) <ex p i 2i z ^ + ^<)> 

N 2 Nl J » VT/ I i=0 J T 

(C2) 

Here F determines the shape of the ramp and satisfies F(Q) - 
and F(l) = 1, and the expectation value (...)t- refers to an 
equilibrium system at temperature r . If we now rescale the 
time and space coordinates in units of A -1 we obtain: 



S 



N a 2N 



(OxPo)"A 
2 N N\ 



t ( 
^' \At) ^ 6XP 2 ' Z e '^' + g ''^ 



V i=0 



> , 

(C3) 



where we define = 0(Ax,, Af,). Due to the scale invariance 
of the underlying theory (and provided that Tq <k //) we have 
that 



(C4) 



validity of the generalized scaling relation (CI 1. This scal- 
ing relation is expected to be valid in the mean-field model 
as well, corresponding to the K — > 00 limit of the one- 
dimensional case. 



Appendix D: Spatially integrated coherence 

Up to this point we considered only the expectation value 
of the local coherence amplitude cos( V20). As explained in 
the body of the article, in order to distinguish between one 
dimensional and mean field dynamics, it is useful to compute 
the full distribution function of the complex operator 



1 r L/2 

LJ-L 



(Dl) 



l-L/2 



Combining this result with Eq.(C3i one can easily show the 



For comparison, in Fig j7] we present the evolution of the full- 
distribution function of Z for (a) the mean-field and (b) the 
one-dimensional cases. In each subplot, the first column con- 
tains a colorplots, proportional to the probability of observing 
a specific complex value of the operator Z. The second col- 
umn refers to the real-phase of the complex operator Z. The 
mean-field and one-dimensional cases look identical for short 
times, but differ at longer times. In particular, the side-peaks 
at (p » +7T disappear as a consequence of the one dimensional 
dynamics. The third column represents the probability distri- 
bution of the absolute value of Z. For the mean-field case, 
these simpliy correspond to delta-function peaks at the maxi- 
mal value of |Z| = 1 . In contrast, for the one-dimension case 
P(\Z\ 2 ) oscillates between the initial Gumbel distribution and 
narrower peaks, indicating the formation of long-range order 
in the relative phase. 
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FIG. 7: [Color Online]Full distribution function of the spatially integrated coherence \Dl\ , according to: (a) Analytic solution of the simple 
pendulum with stochastic initial conditions; (b) Numerical solution of the semi-classical equations in one dimension Both subplots 

represent in the first column a colorplot of the distribution function of Z on the imaginary plane. In the second column, the distribution 
function of the phase of Z. In the third column, the distribution function of the absolute value square of Z. Numerical values- (a) = 0.1, 
time is given in units of the inverse plasma frequency A = 2(j J _/z)~ 1 ' . (b) j±/n = 0.1, K = 25, L = 400, corresponding to the scaling limit 

j±ln 1> K » 1, AL » 1. Time is given in units of the inverse gap A « ((4n/K)j ± po) 



